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' We consider the electron density at the position of an ion of charge Ze in 

' a plasma under conditions approximating those in the core of the sun. Nu- 

merical calculations have shown that the plasma effects on the density, over 

o 

\ and above the ordinary Coulomb factors that are obtained in the absence 

' of electron-electron interactions, are well represented by a reduction factor, 

' exp{— Ze^/9K£)}, where (3 is the inverse temperature and is the Debye wave 

d ' length. Although this factor is the direct analogue of the Salpeter enhance- 

ment factor for the fusion rates in stars, the elementary considerations that 
establish it in the fusion case are not applicable to the determination of the 
electron density and the resulting electron capture rates. We show analyti- 
cally, through a sum rule that leads to a well-defined perturbative approach, 
that in the limit of Boltzmann statistics the Salpeter factor indeed provides 
the leading correction. We estimate residual effects, both from Fermi statistics 
and from short range terms. 
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I. INTRODUCTION AND SUMMARY 



The purpose of this paper is to provide a clear and rigorous treatment of the process of 
nuclear electron capture in a plasma within the mean field approximation. 

It is illuminating to begin by reviewing the effects of the surrounding plasma on the 
fusion of positively charged ions. The experience in looking at these plasma effects under 
solar (or weak screening) conditions is that the correction calculated by Salpeter (1954) is 
considerably larger than any other effect. The essence of this correction is replacement of 
the Coulomb potential by the Debye screened potential. Since the classical turning point for 
the ionic barrier penetration is at a much smaller radius than the Debye radius, it suffices 
to evaluate a screening energy correction to the barrier penetration problem, an energy 
determined by taking the difference between the unscreened and screened potentials as r 
approaches zero. This energy difference is given by /S.E = —e^Z^Z^i^D where Z^ and Z^ are 
the charges of the two nuclei. Then the obvious statistical argument gives the enhancement 
factor exp{/5e^ZiZ2Ki3}. 

However it is not at all clear that this factor applies, as a matter of principle, to the 
plasma corrections to the electron density at a nucleus. One difference is that we now 
deal with an attraction rather than a repulsion, so that at least the language of the above 
qualitative description must be changed. But the important difference is that the thermal 
wavelength for the electron is much greater than that for the ion; for the solar problem it 
is only a little smaller than the Debye length. The WKB approximation gives a theoretical 
justification for considering the outer regions of the electronic cloud surrounding the nucleus 
to be governed by a classical statistical distribution determined by the temperature, chemical 
potential, and local electrostatic potential. But for an electron at a distance from the nucleus 
of one wavelength, the WKB approach is not applicable. Thus the simplest argument for 
the Salpeter factor (now a suppression) seems not to apply to the electron case. 

Nevertheless, in the case of greatest current interest in solar processes, electron capture 
in ^Se, more detailed calculations of electronic wave functions have given screening related 
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reductions to the r — electron density that are only shghtly less than would have resulted 
from a Salpeter formula. In these calculations, beginning with those of Iben, Kalata and 
Schwarz (1967), then sharpened by Bahcall and Moeller (1969), and by Johnson, Kolbe, 
Koonin, and Langanke (1992), the screening induced reduction comes about in an way that 
we describe below. 

1) The Saha equation is used to determine the degree of occupation of the bound states 
in ^Be in the medium assuming pure Coulomb electron-ion forces; for example, there turns 
out to be roughly a 20% probabihty of occupation for both the m = ±1/2 IS states (Iben, 
Kalata, and Schwartz 1967). Leaving out screening, the IS states give a contribution of 
about 35% of the continuum contribution to the r — electron density at the nucleus. 
The higher bound states in this picture have appreciable occupation as well, and contribute 
another 6% of the continuum value. 

2) When Debye screening is introduced, as shown by Bahcall and Moeller (1967), the 
contribution of the continuum states to the r — density is changed by a very small amount, 
of the order of 1%. 

3) The screening reduces the occupation factor for the IS state by a bit; it changes the 
wave function at r = of this state by quite a bit more, with the end result that in the 
screened problem the IS contribution is reduced to about 20% of the continuum contribution. 
As pointed out by Gruzinov and Bahcall (1997), combining this reduction with the complete 
removal of the higher bound states gives a net reduction due to screening that is quite close 
to that predicted by a Salpeter formula. 

Gruzinov and Bahcall address the problem through a quantum diffusion equation eluci- 
dated by Feynman (1990) and applicable in the case of Boltzmann statistics. The approach 
provides both a qualitative understanding of why the Salpeter factor is a good approxima- 
tion in the high temperature limit and a computational framework for incorporating the 
effects of classical plasma fluctuations. Unlike the other approaches described above, it does 
not need to calculate occupation probabilities from the Saha equation. (The Saha equation, 
and indeed the "percentage occupancy" that it describes, are not well defined once electrons 
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interact with themselves.) 

In the present work, which we view as complementary to that of Gruzinov and Bahcall, we 
show analytically that in the Boltzmann limit the Salpeter correction is indeed the leading 
correction in a well defined perturbative approach. Furthermore, in the leading order of 
approximation, Fermi statistics can be maintained at little cost, since the generalization 
of the Salpeter multiplicative factor, in going from pure Coulomb to screened, is a simple 
displacement of the electron chemical potential fie- Thus the Salpeter factor is regained in 
the Maxwell-Boltzmann limit in which everything is proportional to exp{/?/ie}. 

In establishing these results we do not need to separate bound from continuum parts, or 
use a Saha equation. We have also calculated leading small corrections to this basic result. 
In contrast to the fiuctuation corrections of Gruzinov and Bahcall (1997), the additional 
terms depend on the quantum mechanics of the plasma. 

We turn now to describe our results; the mathematical details that support them are 
relegated to the Appendix. 

As in the previous work, we exploit the large nucleus - electron mass ratio so that the 
nucleus may be treated as a fixed point of charge Ze, with the capture rate proportional to 
the electron density at this point which we take to be the coordinate origin. In many-body 
field theory language, this density is given by the thermal ensemble average 

D{p,f,e) = {^p\o)m)p- (1-1) 

HereQ /3 = 1/T is the inverse temperature and /ig is the chemical potential of the electrons. 
This general formula includes all possible corrections resulting from the electron-plasma in- 
teractions to the electron density at the nucleus. The thermal expectation value of the field 
operators is the coincident point limit of the two-point, single-particle electron Green's func- 
tion, including all plasma interactions as well as having an additional fixed point charge Ze 



^We employ energy units for the temperature T and use natural units in which Planck's constant 

n = i. 
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at the coordinate origin. This single-particle Green's function has a standard representation 
in terms of the single-particle irreducible electron self-energy function S. The self-energy 
function may be divided into two parts: 1. A part that contains all the terms depicted 
by all the graphs which end with a single Coulomb photon line connected to the electron 
line. These graphs define an effective, external, local screened potential Vs{r) in which the 
electron propagates. In the small Z, dilute plasma limit, this potential is the Coulomb po- 
tential produced by the heavy nucleus as modified by the plasma polarization accounted for 
by the familiar ring graph sum. In general, however, Vs{r) contain all possible interactions 
of the heavy nucleus with the plasma and all possible plasma interactions, except that all 
these interactions are communicated in the end to the electron by a single Coulomb pho- 
ton exchange. 2. A part that contains all other plasma effects. These entail at least two 
Coulomb photon lines attached to the electron line. The leading correction in this part is 
the exchange energy correction that is familiar in the Hartree-Fock description of atoms. 

In this paper, we shall investigate only the corrections to the electron density at the 
nucleus resulting from the screened potential^] Vs{r). Thus the dynamics of the electron 
field operator ip{r, r) is governed by the screened potential Vs{r) with no other particle- 
particle interactions. The electron field operator satisfies the simple Schrodinger equation 
for a particle in the potential Vs{r). The Green's functions of the theory are of the same form 
as those in the completely non-interacting theory except that in their spectral representation, 
free-particle wave functions are replaced by the corresponding Schrodinger eigenfunctions in 
the potential Vs{r). Hence, with no further approximation, the density reads 



^Although we have given a precise definition of Vs{r) in terms of the single-photon line reducible 
contribution to the electron self-energy function S, the following considerations apply to any local 
mean-field potential Va/ (r). However, any terms that are added to our definition of Vs{r) must 
then be subtracted from the remainder of S, and the net effect of the resulting S must be shown 
to be small. 
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Ds{/3,i^e) = J dEn{E,i2e) \ME;0)f , (1.2) 

where the integration imphes in addition a summation over possible discrete bound states, 
and where 

2 

n(E,/Xe) = ^^(^_^^)^^ , (1.3) 

with the 2 in the numerator accounting for the 2 spin states; ips{E] 0) is the properly 
normalized Schrodinger wave function for energy E. It should be emphasized that the result 
( |1.2|) automatically accounts for the proper weighting of the bound state contributions, no 
additional "Saha-like" reasoning need be done. 

In the limit in which the screened potential is replaced by the Coulomb potential 

Vcir) = , (1.4) 

r 

we have the continuum wave functions for < E < oo with 

in which m is the electron mass, E = p^/2m, = l/e^m is the electron Bohr radius, and 
rj = Z/aop is the usual Coulomb parameter. In addition, there is the infinite set of bound 
state wave functions giving 

\MO)f = ^^, (1.6) 



with the bound state energies 

K = -f^. (1.7) 

We shall denote the the electron density at the nucleus in this Coulomb limit by Dc{P, /ie), 
with 

^ roc 

DciP,fXe) = Y.\MO)fri{En,iXe)+ dE n{E , ^i^) IME] 0)\\ (1.8) 

n=l 
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We shall express the corrections in terms of the electron density (ne) p in the plasma far 
away from the capturing nucleus. Neglecting interacting plasma effects, this density is given 
by 



K)/3 = / j^n{E,fie 



2V2 3^3 



where, in the second line Ae is the electron thermal wavelength defined by 



K = f-^, (1.10) 

V m 

and we have expanded the denominator in the Fermi-Dirac distribution n{E, fi^), performed 
the momentum integrals, and kept the first two corrections to the classical statistics limit. 
For our numerical corrections, we shall use the parameters stated by Gruzinov and Bahcall 
(1997) which describe the solar interior at a distance 6% away from the sun's center. We 
shall also write the parameters in essentially atomic units, except that we shall display the 
units. Thus we take 

p = 0.0215 (ao/e^) , (1.11) 
corresponding to a temperature T = (e^/ao)/0.0215 = 1.27 KeV. This temperature gives 

Ae = 0.368 ao. (1.12) 

The electron density is taken to be 

{n,)f, = 9.10 /al, (1.13) 

which gives by Eq. (|1.9| ) 

^ 0.245 . (1.14) 
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The resulting Debye wave number defined^ by = 47re^/32(ne)/3 reads 

KB = 2.22/ao, k;^^ = 0.451 ao . (1.15) 

We shall first demonstrate that the electron density D£){(3,fie) at the nucleus for the 
Debye potential 

Voir) = —e-^''" (1.16) 

may be expressed, to a good approximation, in terms of the Coulomb limit and then describe 
the correction that arise from a more accurate treatment of the screened potential. This 
relationship between the Debye and Coulomb density follows from a sum rule (|A25|) proved 
in the Appendix^, 

J dE [\i,n{EM'~\ME'-Qf} = if^' (1-17) 

where 

E' = E - Ze^D- (1-18) 
We make use of this sum rule to write 

DDH3,l^e) = J dEn{E,^i,) \ME';0)\^ + Rd{(3, ^le) , (1.19) 

in which 

RD{P,f^e) = J dE[niE,fXe)-niO,iXe)] {\ME;0)\'' - \ME';0)\'} + ^-n{0,iXe) (1-20) 



^Accounting for Fermi-Dirac statistics, the electron contribution to the Debye wave number is 
given by Kj^ g = 4:7:e'^d{ne)i3/dfie- This correct definition reduces the total Debye wave number by 
2%, but it entails only a negligible 1/2% effect for ^Be capture. 

^Note that the convergence of the integration for — > oo is deUcate: The two squared wave 
functions must be subtracted at the same energy; their separate integrals do not exist. 
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will be shown to be a quite small correction. Translating the energy in the integration in 
Eq. ( |1.19| ) to remove the displacement shown in the definition ( |1.18| ) of E' and noting that 



this translation just shifts the value of the chemical potential /Xe in the weight (|1.3| ), we see 
that Eq. ( |1.19| ) may be expressed as 



Dr>{f3, fie) = DciP, fie - Ze^Ko) + Rd{(3, f^e 



:i.2ii 



To keep simple analytic forms, we shall express the corrections in terms of the bulk 
electron density (^e)/3 in the plasma. Including the first non-classical correction, 

1 



-^n(0,/ie) = (ne)/3- 



167ran 



167ran 



1 



1 



2^2, 



:i.22) 



Since the sum rule implies that on the average the difference between \il)s{E;Qi)'f and 
\iIjc{E'] 0)1^ is small, the integral in RD{f3,fie), 

lD{l3fie) = j dE [\^D{E-{))f -\^c{E'-{))\^][n{E,fie)-n{{),fie)] , (1.23) 

should give a small contribution. Moreover, the sum rule has been used to subtract n{0,fie) 
from n{E, fie), and this difference becomes large only at energies that are large on the atomic 
scale, large energies that are on the order of the temperature T. At large energies, the terms 
in the perturbative development of the wave functions converges rapidly and so this high- 
energy contribution can be computed analytically. Said another way, this contribution gives 
the leading term for small /3 and, in view of Eq. (|1 . 1 1|) , P is indeed small in the relevant atomic 
units. For moderate energies, the integration in Eq. ( |1.23|) gives terms that are linear in (3. 
The high-energy contribution obtained in Eq.( A3l|) of the Appendix is of order ^/J3 which 
is larger for small f3. Since the Fermi-Dirac corrections to the Maxwell-Boltzmann limit are 
small, we include only the first correction in this high-energy contribution. Combining the 
result of Eq. ( A31 ) with Eq. (|1.22|) , we obtain the leading correction to the remainder 

1 



J-, \ / \ Zkd^I J (^O^^D 

RD{P,f^e) = {ne),^^-^ 



'3Z aoKD \ / 1 I3k% 
'1 12 J ^ IT 2m 



3V2 



:i.24) 



For the parameters listed above, 

R^(p^ = 0.018 Z{0.28 [1 - 0.16] + 0.098 {Z - 0.25)[1 - 0.26]}{ne)p . (1.25) 

For Z = 4, this gives Rd{P, /^e) = 0.04 (rig)/?, while as we shall see shortly, the corresponding 
capture density R{Pfie) is about 4 {ne)^. Thus this leading additional correction to the basic 
correction provided by the Salpeter factor is only 1%, and any further corrections should be 
smaller yet. 

The Debye potential (|1.16|) that we have been using has the linear term —Ze^Kor in its 
short distance limit. As the work in the Appendix shows, it is this term that gives the non- 
vanishing value to the right-hand side of the sum rule ( P..17D . No linear term in r appears 



in the correct screened potential Vs{r) since such a term would give rise to an unphysical 
screening charge density [ — V^r = —l/r\. However, the corrections that remove this linear 
term at extremely short distances arise from wave numbers in the Fourier transform of the 
potential, Vs{k), that are of order l/A^, where is the thermal wave number of a species s 
particle in the plasma. For the ions in the plasma, the distance is very small in comparison 
with the other relevant distances in our problem, and this cut off in not important. Indeed, 
the work of the Appendix shows that this effect for the ions in the plasma is of order m/M^ 
relative to the small corrections that we have already displayed. But for the electrons in the 
plasma, the cut off at Ae is as important as the other corrections that we have displayed. 
This effect of the electrons in the plasma is computed in the Appendix in the dilute plasma 
limit and to leading order in the small parameter f3Z'^e^ /2aQ. In these limits, the correction 
given by Eq. ( [A41|) reads 

AD(/?,/ie) = (ne)f^/, (1.26) 
87r'^ao 

where / is a parameter integral that has the numerical value / = 1.105. For the parameter 
values used before, this gives a negligible effect: 

AZ}(/3,/Xe) = 0.001 Z(ne);3. (1.27) 
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We have now shown that to within an accuracy on the order of 1%, the screened electron 
density Ds{l3,fie) is given by the simple Coulomb density, but with a translated chemical 
potential, Dc{l3,fie — Ze^no)- It remains to evaluate this main term using Eq. ( |1.8| ). Fol- 
lowing Gruzinov and Bahcall, we define an enhancement factor which is the ratio of the 
electron density at the nucleus calculated in various schemes to the average electron density 
given in Eq. ([L9|) , 

w = ^^. (1.28) 

We denote by ws the complete result of the shifted chemical potential Coulomb density 
including the effect of Fermi-Dirac statistics and including the small correction shown in 
Eq. (|1.25| ). The same quantity, but using Maxwell-Boltzmann statistics, is labeled ws^b- 



The ratio with no screening corrections, the simple Coulomb result is written as wc for the 
case of the full Fermi-Dirac statistics and wc,b in the Maxwell-Boltzmann limit. Finally, 
the results of Gruzinov and Bahcall (1997) will be denoted as wgb- As discussed previously, 
we use the parameters given by Gruzinov and Bahcall in order to compare our results with 
theirs. The ratios given by the various schemes as well as a direct comparison with their 
results are displayed in the following table: 



z 


1 


2 


3 


4 


5 


6 


Ws 


1.38 


1.89 


2.63 


3.67 


5.12 


7.19 


WS,B 


1.39 


1.94 


2.72 


3.84 


5.46 


7.82 


Wc 


1.43 


2.05 


2.96 


4.29 


6.22 


9.04 


WC,B 


1.45 


2.11 


3.10 


4.59 


6.84 


10.3 


Wgb 


1.39 


1.94 


2.73 


3.85 


5.50 


7.90 


WGB 


1.00 


1.00 


1.00 


1.00 


0.99 


0.99 


WS 
WS,B 


0.99 


0.97 


0.97 


0.96 


0.94 


0.92 



The ratio of ratios ws,b/wc,b is essentially the Salpeter factor exp{— /^Zc^kd} which 
varies from 0.95 for Z = 1 through 0.83 for Z = A to 0.75 for Z = Q. The small term 
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RniP, /Ue) displayed in Eq. (|1.25|) gives a correction to ws,b that varies from 0.4 % for Z = 1 
through 1 % for both Z = 4 and Z = 6. The next-to-last row in the table, the ratio of 
ratios ws,b/wgb should be precisely unity if our approximations were without error and the 
computation of Gruzinov and Bahcall were precise.^ The last row in the table, the ratio of 
ratios ws/ws,b, displays the effect of Fermi statistics on the capture rate. 

We conclude that the numbers from the calculations cited earlier are sufficiently accurate 
to determine the electron capture in rate in ^Be to the precision needed for analysis of 
future solar neutrino results. The development presented in the present work first provides 
a rigorous basic formulation that unambiguously describes the electron density at the nucleus 
in the screened field approximation with no need of any considerations of the Saha type. It 
then gives an analytical way of understanding the main features of previous calculations, as 
well as an approach that may be useful in addressing related problems in the future. 
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APPENDIX A: SUM RULE, HIGH-ENERGY BEHAVIOR 



Here we shall derive the sum rule and other results used in the text. This will be done 
by examining the high-energy behavior of S-wave, radial Green's functions defined by the 
inhomogeneous differential equation 

1 £ 



2m dr^ 



Q{E; r, r') = 6{r — r') 



(Al) 



together with outgoing wave boundary conditions. In particular, we shall first compare the 
Green's function Q£){E;r,r') for the Debye potential 



VD{r) 



-Kor 



with the Green's function Qc{E'',r,r') at the translated energy 

,/2 



p_ 
2m 



E' = ^ = E- Ze^KD , 



for the Coulomb potential 



Vc{r) 



This comparison is most easily done by noting that 



/•oo 

gD{E-r,r') = gc{E'-r,r')~ / df gd^E' ■r,f)l^V{f)gD{E-f,r') 

Jo 



where 



(A2) 



(A3) 



(A4) 



(A5) 



AV(r) = Voir) - Vc{r) + E' - E 
■1 



-Ze^ 



(A6) 



This integral equation defines a Green's function ^^(-E; r, r') which obeys the proper differ- 
ential equation ([Al| ) with the Debye potential Voir), and this Green's function is defined 
with outgoing wave boundary conditions if these boundary conditions are obeyed by the 
Coulomb Green's function gc{E';r,r'). 
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The high-energy behavior of the Green's function may be obtained by iterating the 
integral equation ( [A5| ) to form the perturbative series 



poo 

gD{E;r,r')=gc{E';r,r')- dfgc{E';r,f)AV{f)gc{E';f,r') + ---. (A7) 

^0 

The high-energy hmit probes the short-distance hmit of the perturbation, AV{f) — >■ 
—Ze^K,j^f/2, and the high-energy hmit of the Coulomb Green's function which (by sim- 
ple dimensional reasons) has an overall factor ofm/p. The Fourier transform involved in the 
f-integration of the leading f term in AV{f) gives rise, in the high-energy limit to another 
factor of Thus the successive terms in the perturbative development are smaller by a 

factor of in the high-energy limit, and, as we shall see, it suffices for our purposes to 
retain only the first correction as shown in Eq. CKl). 

The Coulomb Green's function may be constructed in terms of solutions to the S-wave 
Coulomb radial Schrodinger equation, the homogeneous counterpart of the Green's function 
differential equation (Kl). These solutions are confluent hypergeometric functions which 



have standard integral representations. The functions we need may be defined by 

roo 

A{E; r) = -2ipr / dt e'P'^^'+^h''''{l + t)'" , (A8) 
Jo 

and 
where 

Ze^m Z / . X 

r/ = = . AlO) 

V P«o 

It is a straight forward matter to check that these integral representations obey the Coulomb 
S-wave radial Schrodinger equation^. It is also not difficult to establish the limits 

r ^ : A{E;r) 1, (All) 



^The application of this differential equations results in integrals of total derivatives that vanish. 
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r ^ : B{E;r) r , (A12) 

and the asymptotic behavior 

r-^oo : A(E- r) -> e'^" ( —] r(l - ir]) . (A13) 

\2pr J 

Thus B{E; r) is the regular solution and A{E; r) has outgoing waves. Moreover, the r — 

limit and the constancy of the Wronskian give 

dA{E:r) , , ,dB(E;r) , 

W = ^-^B(E;r) - A(E;r) = -1 . A14 

ar ar 

Accordingly, the Coulomb Green's function has the construction 

GciE- r, r') = 2m A{E- r>) B{E- r<) , (A15) 

where r> , r< are the greater and lessor of r , r'. 

The results that we need are obtained by examining the high-energy behavior of 

A{E) = hm -[g^(E;r,r')-gc{E';r,r')]^. (A16) 

To the order of accuracy that we need, Eq. ([A7|) and the Coulomb Green's function con- 
struction ([A 151) give 



A(E) = -(2m? drAV(r)A(E;r?. (A17) 

JO 

In the high-energy limit, the Coulomb parameter rj = Z/pao becomes small, and so the 
integral representation (^) gives 

A{E'; r)2 ^ e^'P'"- - Ap'r 7]' H dt 6^'^''^^+'^ In ( . (A18) 

Jo \l + tj 

Moreover, since the high-energy limit involves only the short-distance behavior of the per- 
turbing potential, we may approximate 



1 2 S 2 

— Kr,r Kr,r 

2 ^ 6 ^ . 



(A19) 



To our order, only the first term here contributes in the i]' correction to A(i?';r)^, and we 
have 
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A{E) = Ze\2mf \ + + H In (^) ] . (A20) 

The integral that appears here is reduced to an elementary integral by the variable change 
1 + 1 = 1/ X. Thus, after a little algebra, we find that the leading high-energy limit is given 
by 

in which 



A 



^ (A22) 
2m ^ ' 



12 4 

To obtain the desired sum rule, we note that the Green's functions have a spectral 
representation. Since the radial functions that we have used are related to the total wave 
function by 

V,(^;;.) = !f(^y„o = !i(^_i=, (A23) 

this spectral representation yields 

ACE) ^ r \^s{m?-\^c{E'-M' . , 241 

Air J E-E-ie ' ^ ^ 

where the integration implicitly includes a sum over all bound states. In view of the high- 
energy limit ( |A21| ), we conclude that 

I dE {\ME;0)\'~\ME';0)f} = (A25) 

which is the sum rule used in the text. Since 

Im i . = nSiE - E) , (A26) 
E — E — te 

we also have 

\MEM'-\ME'M'=^^^^, (A27) 
which we shall now use to obtain the corrections presented in the text. 
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The result ( |A27|) expresses the integral in the remainder i?£)(/5,/ie) defined in Eq. (|1.23|) 



as 



dE 



ImA(E) p 
47r2 



(A28) 



Since, in the relevant atomic units, the inverse temperature /3 is quite small, we shall evaluate 
this integral in the small f3 limit, which should give the leading correction. The integration 
over finite energies gives a result that is first order in (3, but, as we shall soon see, the 
integration region of large (3 gives a larger contribution for small (3, a contribution of order 
^/J3. This larger contribution is given by the high-energy limit ([A21j ). Since the departure 
from classical statistics is small, we shall retain only the first-order correction in exp{/?/ie} 
and thus compute 

dE 



Z A 



ZA 



(A29) 



Using Eq. ( |1.9|) to express this in terms of the bulk electron density gives 



with the result ([A22| ) finally yielding 



2nan 



3^2 



IT 



3V2 



3Z aoKD\ / 1 l3Kj-) 

1 i2~; V 



(A30) 



(A31) 



We turn at last to examine the correction to the electron density at the nucleus that comes 
from the difference between the Debye potential VD(r) and a more accurate screened poten- 
tial Vs{r). Since the two potentials have the same long-distance behavior, their difference 
becomes important only at high energies where the effects of the Coulomb interaction be- 
come small. Hence, a good estimate is obtained by replacing the Coulomb function A{E'] r) 
by the plane wave exp{ipr} in Eq. ( [A17] ): 



AAE) = -(2m? / dre^'P' {Vs(r) - VD(r)} 
Jo 



(A32) 
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The Fourier transform of the screened potential may be written in the general form 



A;2 + 47rn(A;) 



(A33) 



With 47rn(A;) = K^j, this gives the Debye potential. Hence 
Vs{r) - Voir) = Ze^ 



(dk)^^,^ 47r 47r [n(/c) - n(0)] 



8Ze' 



°° dk sin kr . 

[u{k) - n(o)] , 



(A34) 



k"^ kr 

where is the second approximate equality we have performed the angular integration and 
neglected the denominator corrections that are of order nj^ and do not contribute to the 
leading high-energy behavior. Using 



/ 

^0 



°° dr 



e^'^'' sin/cr 



In 



2p-k 
2p + k 



(ASS) 



and thus obtain 



ImAi(E) = 



16mZ r"^ dk 



- [U{k) - n(0)] In 



2p-k 



2p + k 



(A36) 



ao Jo k^ 

Again for simplicity using Maxwell- Boltzmann statistics, this gives the correction to the 
electron density at the nucleus of 

p2 I imAi(E) 



AD{P,^^)=2eP^^ r^exp(-/5|^| 
Jo m [ 2m] 



47r2 



(A37) 



We now represent the plasma polarization function Ii{k) by its one-loop, ring graph ap- 
proximation. Using the results of Brown and Sawyer, their Eq's. (A41) and (A42), we find 
that 

13k'' 



Ti{k) - n(o) = E/'/^K)/? /' 

t/ 







r dt 


exp 1 


Jo 





2M. 



■t{l - 1) 



(A38) 



where the sum runs over all species s of the particles in the plasma. We perform the 
integration over p by writing k — 2px to obtain 



dt 



In 



1 — X 
1 + x 



m 



In l + AxHil-t) — 



(A39) 
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The ratio m/Mg is very small for the ions in the plasma. In the ionic case, the final 
logarithm above gives only a term of order m/Mg. Thus only the electrons in the plasma 
give a significant contribution to this high-energy correction. Since the electron density is 
half the total number density in the completely ionized plasma, we may write 

1 



(A40) 



The t integration in Eq. (|A39|) is elementary, and 



AD(/?,/x,) = K)^^/ 



in which / is the analytically intractable integral 

dx 



Jo X* ^ 



x^ In 



X 



x| In 



1 — X 



1 + x 



(A41) 



(A42) 



whose numerical integration gives / = 1.105 
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